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DNA's shape mostly lies on its total energy F. Its corresponding equilibrium shape equations can 
be got by classical variation method: letting the first energy variation S"'F = 0. Here we not only 
provide the first variation S^F but also give the second variation S^ 2 'F in planar case. Moreover, 
the general shape equations of DNA are obtained and a mistake in [Zhang, et al. Phys. Rev. E 70 
051902 (2004)] is pointed out. 
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1. Introduction 

Closed duplex DNA molecules often have complex configurations, such as supercoil, and those structures play an 
\ important role in gene regulation. 1 Theoretical analysis about those configurations are based on the elastic rod theory. 
In 1859, Kirchhoff given the equilibrium shape equations of thin elastic rod and his methods has been generally used to 
describe the holistic conformations of duplex DNA. 2-5 Benham 2 treated a homogeneous isotropic rod of the DNA by 
Kirchhoff analogy and gave some possible shapes. As for the closed duplex DNA molecules, its equilibrium equations 
are obtained 6 , and considerable results are attained. 7-11 

Following the development of modern experimental techniques, such as optical tweezers, 12 micromanipulation, 13-15 
and other techniques, 16-18 researchers are presented with more opportunities to probe into the microstructure of 
individual great molecule than ever before. The ever growing volume of experimental data provides us with the 
I < probabilities to improve the existing theoretical models. 

Recently, considering the microcosmic configurations of DNA, more practical models, such as the anisotropic 
model 19 are put forward, and new energy terms are introduced, 20-24 and new methods, such as Monte Carlo 
simulation 20,25 and path integral, 21 are adopted. Due to those improved theories, some theoretical results are highly 
consistent with the experimental values. 21,24 

The geometrical configuration of duplex DNA is shown in Fig.l. The axis of the molecule can be written as R(s) 
parameterized by the arclength s. The total energy of a simple elastic model for the closed duplex DNA can be 
written as 9 
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where K = K{s) and r = r(s) are the curvature and torsion of the axis R(s). A and B are the bending rigidity and 
torsional rigidity, respectively. Some researchers also choose a simpler model without the torsional term. 26 

This letter is composed as follows. In Sec. 2, we will get the equilibrium shape equations of closed duplex DNA by 
variation theory in its nature coordinates system. In Sec. 3, a general model is considered and its equilibrium shape 
equations are attained. In Sec. 4, we give the second variation S^F in planar case. Finally, a short discussion is 
given in Sec. 5. 



2. Classical Equilibrium Shape Equations of Closed Duplex DNA 

For a more general case, we choose the total energy of DNA 

F = j ^(K -C ) 2 + B(t -T ) ds + A j> ds, (2) 

where Cq and To are two constants, and we call them spontaneous curvature and spontaneous tension respectively, 
and A is a Lagrangian multiplier which can be taken as the tensile rigidity. Fig. 1 shows a fragment of R(s) in the 
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FIG. 1: The geometrical configuration of duplex DNA molecule. The central line of the ribbon follows the molecular axis with 
coordinates R(s) parameterized by the arclength s. A nature coordinates is defined: the tangent vector a, the main normal 
vector n and the vice-normal vector f3 



orthogonal coordinates system (x,y, z). Meantime, a nature (local) coordinates can be described as 

(a,n,/3). 



(3) 



Where, a = R is the tangent vector (an overdot denotes differential with respect to s), n = R/K is the main normal 
vector and (3 = a x n is the vice-normal vector. Between those unit vectors and K and r, there are the Frenet 
formulacs 



a = K(s)n, 

h = -K(s)a + t(s)/3, 

= — r(s)n. 

The equilibrium shape equations will be extracted by analyzing the first variation 

5V>F = 0. 

Letting g = R x ■ R x (R, x = ^), we have ds — y/gdx. So, Eq. (7) induced 



(4) 
(5) 
(0) 



(7) 



A(K - C Q )6 (1) K + B6 {1) t] Jg 



A 



(K - C ) 2 + B(t - T ) + A 



<5 (1) (Vs)}^ = 0. 



(8) 



Before we get the equilibrium shape equations, some useful identical equations should be calculated (see appendix A). 
The R(s) under small perturbations can be written as 



R (s) = R(s)+ip{s)n + (p{s)/3, 
where i/)(s) and ip(s) are two small smooth functions. We get 

5R = R — R = ipn + (p/3 

5R. X = (SR). X = ^/g(6R), s = V5(^n + ifjh + cp/3 + tpfr) 



(9) 
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and 



Sg = (R, x + 5R, X ) ■ (R, x + SR, X ) - 5R, X ■ 5R, X 

= 2R X ■ SR, X + 5R, X ■ 5R, X = g[2R ■ {SR), S + (SR), S ■ {SR), S 



-2gKip + g 
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Consequently, we have 
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-giKiP + ^gi ip 2 + (ip 2 + ip 2 )r 2 + ip 2 + 2(ipip - ip<p)r +0(3), 



Sg 2 = (g + Sg) 2 - g 2 =g 



1 ( Sg\ , 3 (Sg 
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g 2 Kip - -g * 



2 V 9 J 8 V 9 

ip 2 + (ip 2 + <p 2 )r 2 +ip 2 - 2iP 2 K 2 + 2(ipip - ipip)r 



+ 0(3). 

Where 0(3) means the third and higher orders of ip and ip. We can also attain 

SR = S(g^R, x ) = (g=r + Sg^)(R, x + SR, X ) R 
= (l + ghg^)(5R) tS +g^RSg^ 
This conclusion can be generalized to an arbitrate function V = V(s), which satisfies 

SV = (1 + ghg ^)(SV), S + g^VSg^r 
Submitting (11) and (14) into (15), we get 

SR = KipR + ipn + ipii + ipfi + <pj3 + Kip(ipn + iph + <p/3 + ip'p) 

-^R [ip 2 + ip 2 + (i/j 2 + <p 2 )r 2 - 2ip 2 K 2 + 2{ip<p - ^)t] + 0(3). 

Using (16), we have 

5R=(l + gt5g = r)(5R),,+gtR5g = r. 
Insetting (14) and (17) into (18) we get 

SR = 2KipR + R(Kip + Kip) + ipn + iph + 2iph + <p(3 + ip'/3 + 2<p0 



+KiP 
-R 



2ipn + 2iph + 2ipf3 + 2ipf3 + Aiph + Atpfi + KipR + (Kip + Kip) A 



ip 2 + ip 2 + {ip 2 + ip 2 )r 2 - 2K 2 iP 2 + 2(ip(p - ip<p)i 



+ (Kip + Kip) (ipn + iph + ipf3 + ip(3) 



R 



Tf(ip 2 + if 2 ) - 2KKip 2 - 2K 2 ipip + f(ip(p - iptp) + r 2 (ipip + (pip) 



+ipip + <p>Cp + r(ip(p - ipip) + 0(3). 
Noting K 2 = R ■ R, we have 

S^K = (R-S^R)/K. 
Submitting (19) into (20) and using the identical equations in appendix A, we attain 

S (1) K = (K 2 -T 2 )iP + iP-iip- 2rip. 

Now, we calculate 
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Submitting (19) and (21) into (22), we get S^h. Then, using the identic equations in appendix A, the following 
useful term can be attained 



h-S^h = K 3 ip + K \T 2 ip - Tip - rip + ip 



+K- 1 t 



Tip — T 2 lp + 3tlp + 2t(V> — Tip) + if 



(3) 



+K 2 Kt (t 2 : ip - Tip - 2t0 - Ifij 



Considering h • h = K 2 + t 2 , we have 



5™T = 



h- Pn- KS^K 



ip - ip 



Taking (21) and (23) into (24), we get 

S^ T = K{2T%P + tp) + K- 2 k(r 2 tp-T%P-2Ti; 

+R- 1 fip + 'Hip -T 2 p + p {3) + 2r(ip- pij 

Submitting (13), (21) and (25) into (8), we attain 

'A 
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-Kip 



-{K -Co) 2 + B(t-T ) + X 



+A(K-C ) (K 2 - t 2 ) ip - fp - 2tp + ip 
+BK(2Tip + <p)+ BK~ 2 K (r 2 tp -Tip- 2 T ip - <pj 



+BK' 1 



Tip + 3fip — T 2 ip + p^ + 2r (ip — iip^j 



ds. 



If perturbation is only on the main normal direction (p = 0), we have 



-A(K 3 + 2C t 2 - KC 2 - 2Kt 2 ) -K{\-Bt- BT ) 
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+BK- 2 (Kt- Kt), 
BK- 2 (3Kt -2Kt), 
A(K - C ) + 2BK~ 1 t, 



ipds. 



(23) 
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(27) 



(28) 



(29) 



S W F = <j> (Foip + Fnp + F 2 ^j ds. 

For two smooth functions Q = Q(s) and P = P(s), there is 

/*■>«.-(-!)./*■>«*. 

where n = 1, 2, 3 • • • . Using (28), (27) can be changed into 

5&F = <j> ^AK 3 - ^K (AC% + 2\- 2BT - 2Bt + 2At 2 ) + A (c t 2 + if) 

Because ip is an arbitrary function, so the equilibrium shape equation under main normal perturbation is 

AK 3 - K{ACl + 2A - 2BT - 2Bt + 2At 2 ) + 2A(C q t 2 + K) = 0. (30) 

Another equilibrium shape equation under vice-normal perturbation can be obtained through the former way. Choos- 
ing ip = 0, (26) can be simplified as 

F = BK- 2 Kt 2 + 2BK- 1 tt - A(K - C )f, 

F 1 = K{B -2At) + 2AC q t - BK^t 2 , 

F 2 = -BK- 2 K, 

F a = BK~\ 

S {1) F = <j> {Fop + F^ + F 2 (p + F 3 V (3) ) ds. (31) 
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The corresponding equilibrium shape equation obtained by Fq — + 4-1% — ^L^l = is 

A(K - C )t -(B- 2At)K = 0. (32) 
Eq. (32) can be changed into the following state: 

(K-C ) 2 (2r-Q) = C, (33) 

where Q = B/A, and C is a constant. When Co = To = A = 0, similar equations are attained and an elliptic function 
solution is found in Ref. [6]. In particular, if R(s) is a planar line (t = B = T = 0), only Eq. (30) is valid and is 
reduced to an elliptic different equation 

2AK + AK 3 - {ACl + 2\)K = 0. (34) 

A similar equation was discussed in vesicle research. 27 

Here, two special cases should be taken into account. First, there is a ring solution to Eq. (34): 



K = l/r= ^C 2 + 2\/A, (35) 

where r is the the ring's radius. If there is stable planar ring DNAs, we need 

Cl + 2X/A > 0. (36) 

Second, it might have isolated inflexion points in the two shape equations at s = So, if 

K(s Q ) = C . (37) 

When Co = 0, this problem is discussed in Ref. [6]. At the inflexion points, DNA's shape may has uncertain behavior. 
Let's see an example. Choosing C = in (33) and supposing R(s) is a solution of Eqs. (30) and (33), then, one can 
find that if there is a point s = so in R(s) satisfies K(so) = Co, Eq. (33) will lose the constraint of r and there is 
only Eq. (30) which needs to be satisfied. In this case, Eq. (30) is easy to be satisfied, because it has two variables: 
K and r. Therefore, we can expect that DNA's shape changes variably at its inflexion points. 

3. General Shape Equations 

It is known that the curvature K and torsion r are two characteristic parameters which can determine a line's 
shape. Recently, new energy terms are introduced 20-24 in DNA's total energy. Those new energy terms will influence 
DNA's shape through changing its local curvature and tension. Moreover, DNA's total energy can be written as the 
sum of the elastic energy of its two molecular lines: Fj = §[AK 2 (si) + BT 2 (si)\ds (i = 1, 2) (see Ref. [4]). Where Sj 
is the arclength of one of DNA's molecular lines. Clearly, we can choose one of DNA's molecular lines to determine 
its shape and attain its shape equations by letting S^Fi = 0. 

Here, we suppose that those new energy terms can be expressed by curvature and tension, and we choose the total 
energy of closed DNA 

F=If(K,r)ds, (38) 

where f(K,r) is a function with two variables K and r. Similar method is adopted in vesicle research. 28 ' 29 By the 
former way we have discussed the first order variation of total energy, that is 

S^F = j[{f, K S (1) K + f tT &V>T)y/g + f5W(Jg)]dx, (39) 

where f <K = §t. Submitting (13), (21) and (25) into (39), we get 

F - -Kil>, 

Fi = (K 2 -T 2 )ip-i(p- 2rip + i>, 

F 2 = K(2ti\) + if) + K~ 2 k (r 2 ip -rip- 2rip - <p) 

+K~ X 3TTjj-T 2 ip + fij + (p (3) +2T(ip-Tipj , 
S (1) F = f(F f + F 1 f tK + F 2 f tT )ds. (40) 
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Perturbation only on main normal direction (tp = 0) induces 

F = (K 2 -T 2 )f, K + K- 2 (2K 3 T-KT + Kt) f, T -Kf, 

F 1 = K- 2 (3KT-2KT)f. T , 
F 2 = f, K + 2K- 1 rf, T , 

S^F = J) (f q iI> + Flip + F 2 , if S j ds. 



The corresponding equilibrium shape equation is 

d 2 



K- 2 (^Kt-2Kt)1, 



+ (K 2 - t 2 ) f.K + K- 2 (2K 3 t -Kt + Kt) f, T - Kf 
If perturbation is on vice-normal direction [tp = 0), (40) is simplified as 

F - K- 2 T(KT-2KT)f, T -ff, K , 

F 1 = K- 1 (K 2 - t 2 ) f tT - 2rf, K , 

F 2 = -K- 2 K.f, T , 

F 3 = K- 1 ^, 



0. 



S^F = I (F <p + Flip + F 2 ip + F 3 <p^ ds. 



The corresponding equilibrium shape equation is 



^ («-*f.*) + £ (K-»kf. r ) + ± [K-i (K 2 - r 2 ) f, T - 2rf, K ] 

-K- 2 T(KT-2KT)f, T + if, K = 0. 



(41) 



(42) 



(43) 



(44) 



Eqs. (42) and (44) are the general equilibrium shape equations for closed DNA. Comparing Eqs.(41) and (43) with 
Eqs.(6) and (7) in Rcf. [30] respectively, we find that Eq. (42) is equal to Eq. (6) in Ref. [30] (a(s) = a s (s) = 0). 
But Eq. (43) is not equal to Eq. (7) in Ref. [30]. Acutely, the mistakes in Ref. [31] induce the incorrect result in 
Refs. [31] and [30] (see appendix B). Moreover, similar correct equations are also shown in Ref. [32] (Eqs. (16) and 
(17)), in Ref. [33] (Eqs.(l) and (2)) and in Ref. [34] (Eqs.(77) and (78)). 



4. Second Variation in planar case 



We know that stability condition needs the second variation S^F > 0, so we want to get S^F by the former way. 
For simplify, we only consider the planar case with the energy F — §{K 2 + X)ds. Consider K 2 = R R, there is 



SK 2 = 2R <5R + SR SR. 
Submitting (19) into the above equation, we get 

SK 2 = 2R 5R + 8R SR 

= 2K(K 2 t/j - T 2 ip -fip- 2r<p + i>) + 3ifV + r 4 ((fi 2 + 2 ) 
+ (ip 2 + 2 )f 2 + 4r a (ip(p - ipip) + 2kKip(ip - np) 

+2i(ip<p - 4>ip) + ip 2 + ip 2 — K 2 T 2 (<p 2 + 80 2 ) + Giipip 
+ l6Tipip + 2ip 2 -ip 2 - 6-0-0 + 2r 2 (2-0 2 + 20 2 - -00 - (pip) 
+4r f(00 + tpip) + tp/p - rpip +0(3) 



(45) 



(46) 
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Allying (13), we get the second variation of DNA's energy (note r = ip = 0) 



(K 2 + \)5 {2) ^g + ^g8 [2) K 2 + S^^gS^K 2 dx 
-(3K 2 + \)i> 2 + AK 2 ^ + K 4 tfj 2 + 2KKipip + ip 



ds. 



(47) 



The above equation can be used to analysis the stability of DNA in planar case. For instance, it's easy to find that 
straight line K — is a stable solution of Eq. (34) if A > 0. 



5. Discussion 



In Sec. 2, we give the DNA's shape equations which contain two special constants: Cq and To. But what are their 
physical meaning? Maybe we can take them as the influence of the environmental factors or other complications. For 
example, compared with vesicle research, those two constants can be attributed to the donation of electric potential 
on DNA. 

In (47) we only give S^F in planar case, if we want to get 5^F in a more general case, S^K and S^t should be 
obtained formerly. For instance, if we consider 



SX n = nX^SX 



2! 



-(sxy 



(48) 



where X — X(s) is a scalar function of s, and n is a constant, we can attain all SK n terms by using Eq. (46). 



Moreover, using (16), all SK^ {K {r 



d m K 
ds 11 



rn 



1, 2 ■ • ■ ) terms can be obtained. Finally, because attaining 5t^ 



needs a tedious calculation, we will show it and discuss 8^ 2 'F in another paper. 
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APPENDIX A 

In this appendix, we show some useful identical equations which can be attained easily. 



ri 


R = -K, 


(Al) 


ri 


R = 0, 


(A2) 


ri 


R (3) = R 3 + T 2 R ^ 


(A3) 


ri 


R (4) = 3KK 2 + t(2Kt + Kf), 


(A4) 


ri 


P = T, 


(A5) 


ri 


(3 = 0, 


(A6) 


ri 


(3 = -t 3 - K 2 t, 


(A7) 


ri 


/3 (3) = -KKt - 2K 2 t - 3r 2 f, 


(A8) 


ri 


ri = K 2 + r 2 , 


(A9) 


ri 


ri = KK + tt. 


(A10) 


ri 


n (3) = KK + tt- (K 2 + t 2 ) 2 , 


(AH) 


ri 


h = k 2 + r 2 + (K 2 +r 2 ) 2 . 


(A12) 



APPENDIX B 



In this appendix, we point out a mistake in Ref. [31]. Comparing this letter with Ref. [31], we find that they 
are similar work. Here we use ip and (p to determine the perturbations on the main normal direction n and the 
the vice-normal direction (3, respectively. (Actually, there should be perturbation on the tangent direction a, but 
we find that it will not induces any new shape equation, so we ignore it deliberately.) In Ref. [31], there are three 
corresponding perturbation functions: E\ on the tangent direction a, £ 2 on the main normal direction n and £3 on 
the vice- normal direction (3. Clearly, there are the relationships 



ip = e 2 , <p = £3- 



(Bl) 



Comparing Eq. (25) with Eq. (2.27) in Ref. [31], we find they are not equal. There are two mistakes in Eq. (2.27) in 



Ref. [31] that the first term K £3 in the first bracket should be changed as t £3, and the — c 3 
The correct S^r in Eq. (2.27) in Ref. [31] should be 



(3) 



term should be e 



(3) 



5«T = -^jA-(T 2 £ 3 -f£ 2 -2T£ 2 -£ 3 ) 



K 2K z te 2 + f(-2r£ 3 + 3£ 2 ) + (A - t 2 )e 3 + te 2 + 2t£ 2 + e y 3 



(3) 



(B2) 
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One can find that the above equation is equal to Eq. (25) (note E\ = 0). It is this mistake that induces the incorrect 
shape equations: Eq. (2.32) in Rcf. [31] and Eq. (7) in Ref. [30]. 



